NASA Technical Memorandum 89048 


2nd Annual 



PtrossiaO flags off a ^7afffisQistP 
Sotneairaifling aQia ...... 



DESIGn CHHLLEI1GE 

Held at 



NASA Langley Research Center 
Hampton , Virginia 23665 
December 9-J0 r 1985 


NASA 

National Aeronautics and 
Space Administration 


Compiled by Larry Taylor 

(NASA-TM-89048) fBOCEEDIKGS Cf THE 2ND 
ANNUAL SCOLE N0BKSH0P (NASA) 268 p CSCL 22B 


Langley Research Center 

Hampton, Virginia 23665 

*Spacecraft Control 


N87-1C945 


Onclas 

G3/18 43975 

Laboratory Experiment 



The 2nd Annual &C01& Workshop 

at the NASA Langley Research Center, Hampton, Virginia 
Building 1 1 92 C, Room 1 24 


Monday, December 9 

8:30 Welcoming Remarks John Dibattista 

8:45 Definition of the SCOLE Design Challenge Larry Taylor 


9:00 Panel on Issues in Modeling the SCOLE Configuration 

A.V. Balakrishnan 
Peter Bainum 
Jer-nan Juang 
Leonard Meirovitch 
Yogendra Kakad 
Tony Hotz 
Larry Taylor 

12:00 Lunch 

1:00 Guest Investigators Role in the COFS Program Tony Fontana 

1:15 Panel on Issues in Controlling the SCOLE Configuration 

A.V. Balakrishnan 
Peter Bainum 
Shalom Fisher 
Dan Minnick 
Gene Lin 
Sahjendra Singh 
Steve Yurkovich 
Mike Barret 

5:00 Attitude Adjustment 
6:00 Supper Without Speeches 


Tuesday, December 10 

8:30 Laboratory Tour in Building 1232 

SCOLE Experiment Demonstration .Jeff Williams 

Grid Experiment Demonstration Ray Montgomery 

Optical Sensing Demonstration Bob Bullock 


12:00 Lunch 


1:00 Laboratory Tour in Building 1293B 

Flexible Beam Experiment Demonstration Lucas Horta 

15 Meter Hoop/Column Antenna Structural Tests Keith Belvin 

2:30 End of Workshop 



Table of Contents 

Page 

The SCOLE Design Challenge - An Overview 1 

Larry Taylor 

SCOLE Equations of Motion - A New Formulation 1 4 

Suresh Joshi 

Dynamics of Flexible Spacecraft (SCOLE) During Slewing Maneuvers 26 

Yogendra Kakad 

Functional Analysis of the SCOLE Problem 38 

A. V. Balakrishnan 

Maneuvering of Flexible Spacecraft 50 

Leonard Meirovitch 

Model Reference Control of the SCOLE 62 

Dan Minnick and Howard Kaufman 

Rapid Pointing and Vibration Control of the SCOLE Configuration 82 

Jiguan Gene Lin 

A Perspective on the Control of Flexible Spacecraft 86 

Mike Barrett and K. W. Lips 

Equations du Mouvement d'une poutre Flexible en Rotation Autour d’un Axe 
Lionel R. Passeron 102 

Flexible Beam Simulation 1 13 

Shalom Fisher and Tom Posbergh 

Issues in Modeling and Controlling the SCOLE Configuration 132 

Peter Bainum, A.S.S.R. Reddy and Cheick Modibo Diarra 

Decentralized Control Experiments on a Flexible Grid 217 

U. Ozguner and S. Yurkovich 

A Mathematical Problem and a Spacecraft Control Laboratory Experiment 

(SCOLE) Used to Evaluate Control Laws for Flexible Spacecraft (Reprint) 

L. W. Taylor, Jr. and A. V. Balakrishnan 236 


List of Attendees 


263 



The SCOLE 
Design Challenge 

by 

Larry Taylor 
NASA Langley 
Research Center 
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SPACECRAFT CONTROL LAB EXPERIMENT 
(SCOLE) 
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Benefits of the 

$C01C Program 

O Comparison of Design Approaches 

O Novel Control Techniques 

O Validation of Design Methods 

O Precursor to Co- Investigators for 
Flight Experiments 
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Objectives of the 



0 To Validate Methodologies for the Control 
of Flexible Antenna Structures 

0 Identify Advantages and Disadvantages of 
Various Systems Identification Techniques 

0 Identify Advantages and Disadvantages of 
Various Control Synthesis Techniques 

0 Identify Outstanding Identification and 
Control Candidates for the COFS II Flight 
Experiments 
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flexible 

BE/im-l'ke 

TRUSS 



Justification for tho 

Program 


0 Methodologies for Control of Flexible Spacecraft 

Have Not Yet Been Validated 

* 

0 Methodologies for On-Orbit Modelling of Structural 
Dynamics of Flexible Spacecraft Have Not Been 

Validated 

0 The High Cost of On-Orbit Testing Needed to 
Support Initial Validation May be Alleviated 
by Using Hi -Fidelity Ground Based Experiments 

0 A Direct Comparison of Different Design Approaches 
for the Control of Flexible Spacecraft is Needed 

0 On-Orbit Systems Identification Techniques 
Need to be Compared in Preparation for Actual 
Flight Application 
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Timeliness of the 

Program 


0 Ground Demonstrations on Two Dimensional 
Structures are Nearing Completion - A 3- 
Dimensional Problem is Needed to Advance 
Our Active Structural Control Capability 

0 The Planned COFS Flight Tests Demand 
Development of Our Most Sophisticated 
Flexible Structure Control Capability - 
This Can Only be Accomplished by Ground 
Testing on Structures with 3-Dimensional 
Characteristics and Large Numbers of 
Inertial Sensors and Actuators 
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Results of the 


Program 


O 

O 

O 

O 


NASA- IEEE Design Challenge 

June '83 

Eight Grants/ Contracts Underway 

Feb '84 


SCOLE Workshop 

Dec ‘84 

Experimental Apparatus Ready 

July '85 
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Milestones of the 

Program 
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SCOLE 

Equations of Motion 
A New Formulation 

by 

Suresh Joshi 
NASA Langley 
Research Center 


SCOLE equations of motion are derived 
in the Shuttle-fixed coordinate system 
using a Newtonian formulation. 
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Dynamics of Flexible 
Spacecraft (SCOLE) 
During Slewing 
Maneuvers 

by 

Yogendra P. Kakad 
University of N.C - Charlotte 

Charlotte, N. Carolina 
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The velocity in the body frame is 


v co 







The velocity of the point of attachment in the body frame is 


Vo =. V +• A 


where r is the vector from mass center to the point of 
attachemnt . 


If the body-fixed frame is located at the point of attachemnt, 
a is the position vector of a mass element on the beam from 
the point of origin before deformation. 


The displacement vector is 






O 


The position vector after deflection is 


Tl * 4 - cL. 
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The kinetic energy in the beam is 
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Equations of Motion 


Lagrange's Equations 
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Translational Equation 
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Vibration Equation 
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Functional Analysis 
of the SCOLE 
Problem 

by 

A. V. Balakrishnan 
U. C. L. A. 

Los Angeles, California 
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Figure 8.5 Mode 3, FACT = 0.25, 
0.033353 Hz 
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Table 8.4 Natural Frequencies (Hz) 



Equi librium 

Position 

i 

Zero Gravity 

(Space) 



Massless 


Rigid 

Mode 

Actual 

Cable 

Actual 

Antenna 

1 

0.0 

0.0 

0.0 

0.0 

2 

0.02661800 

0.02661800 

0.0 

0.0 

3 

0.03335345 

0.03335345 

0.0 

0.0 

4 

0.27527161 

0.27550803 

0.0 

0.0 

5 

0.27528201 

0.27551848 

0.0 

0.0 

6 

1.02135153 

1.02136222 

0.95626111 

0.95781737 

7 

1.09274338 

1.09275532 

1.02205468 

1.02304393 

8 

2.86993771 

2.86993909 

2.85798288 

2.91078995 

9 

4.14099978 

4.14100682 

4.12238565 

4.77249201 

10 

7.21202381 

7.21203514 

7.13573328 

7.53400307 

11 

11.8888566 

11.8892693 

11.8067296 

14.3285908 

12 

11.9975825 




13 

11.9979985 




14 

14.5271386 

14.5271319 

14.4703039 

18.0655372 

15 

17.8807340 

17.8965361 

0.0 

0.0 

16 

23.9951935 




17 

23.9951971 




18 

29.4618566 

29.4618689 

29.3765971 


19 

31.9790379 

31.9790406 

31.8650183 


20 

35.6083428 

35.6083908 

35.5681068 



Table 8.5 

Cantilever 1 

Natural Frequencies (Hz) 



Roll Angle (Deg) 

Mode 

0 

180 

1 

0.90329379 

0.74059481 

2 

0.92231250 

0.74130235 

3 

2.85838001 

2.84157907 

4 

4.11631152 

4.06401370 

5 

7.18829155 

7.03295932 

6 

11.8569533 

11.6900310 

7 

14.5078152 

14.3982264 

8 

29.4282394 

29.2763643 

9 

31.9410827 

31.7363829 
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MODEL REFERENCE CONTROL OF LUMPED LINEAR SYSTEMS 
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SPECIAL CASE (PMF) 
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MODEL REFERENCE CONTROL OF DISTRIBUTED SYSTEMS 



II 



CONTROL PROBLEM FORMULATION 



11 


B<3C ) e 







ILLUSTRATIVE APPLICATION TO SCOLE ROLL 



IF 


FOR ILLUSTRATIVE PURPOSES WE WILL CONTROL 



76 


a<5(s- 


REFERENCE MODEL 


! 


i 



CO 

II 


• cr 


cr 

*— 

o 

ii 



cr 

CM 

O 

II 


CM 



cr 

co 

o 

II 

CO 

>f 


I 


77 




RECALL THAT FOR TRUE STABILITY WE NEED 



79 


CURRENT EFFORTS 



80 



PLANNED ACTIVITIES 



81 


Rapid Pointing and 
Vibration Control 
of the SCOLE 
Configuration 

by 

Jiguan Gene Lin 
Control Research Corp. 
Lexington, MA 


82 


STRUCTURAL VIBRATIONS 


o 10 -mode ABC-model (Dr, Joshi) 

o Linearized Dynamics 
-- Decoupled with RB motions 

— Excitation by the BPB LOS Poinging Slew 

O^SMALL' VI BRAT IONS 

— OUTPUT AT BOTH ENDS OF MAST 

o Need Output from the 12 sensors along the mast 

—ACTUAL VIBRATIONS OF THE MAST 

o Need Complete dynamics (RB & EB) 

—Coriolis coupling could further excite Vibrations 
due to resulting moment spillover. 



RAPID LOS POINTING 


o Initial Alignment and BPB pointing slew 

o Coriolis coupling (w X) 

— Moving Shuttle-body-fixed coordinates 

o Direct torque/force spillover 
—Single Axis BB or BPB 

= ^1 ^DEFLECT 

T T T 

I'E" T4 “ ^DEFLECT '1 



VIBRATION CONTROL 


o Control Influences 
-- Actuators at ends of mast only 

— High feedback gains 

— Hence, high additional torques and forces 


o Need Actuation along the mast, 
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• Computational Complexity 
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SYSTEMS IDENTIFICATION EXPERIMENTS (CONTINUED) 
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#Assess adequacy of MLE identification methodology for providing 
reliable structural models for control design 
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Control Problem Issues 
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Robust Control Experiments 
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• Correlate actual results with analyses 



Robust Control Experiments (Continued! 
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•Assess adequacy of control design methodology for 
meeting control objectives 


Equations du Mouvement 
d une poutre Flexible 
en Rotation Autour 
d'un Axe 

by 
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L. FA SSERoN 


EQUATIONS DU MOUVEMENT D'UNE POUTRE 
FLEXIBLE EN ROTATION AUTOUR D'UN AXE 


On determine, dans ce qui suit, les Equations du mouvement 
d'une poutre flexible encastrAe A sa base sur un axe verti- 
cal en rotation. 

On neglige la flexion verticale pour ne s'intAresser qu'aux 
raouvements de traction-compression et de flexion horizonta- 
le j autrement dit, la poutre est supposAe se d§placer dans 
un plan horizontal. 


2.1 - NOTATIONS ET HYPOTHESES 

2.1.1 - Notations 

On denote respectivement : 

- par I = (0,f,5,&) un repdre inertiel, le vecteur £ Atant 
dirigA suivant la verticale, 

- par R = (0,T,^,l£) un repAre liA A la barre considArAe 
comme rigid©, 

- par 0 1' angle de la rotation autour de K permettant de 
passer du repAre I au repAre R. 

La ligne moyenne de la poutre A l'Aquilibre est dirigAe 
suivant l'axe i. 


On dAsigne respectivement : 

/ » 

- par G (x) le point courant de la ligne moyenne, 

- par p(x) la masse volumique de la poutre 

- par S(x) la surface de la section droite de la poutre 

- par I z (x) le moment quadratique de S(x) par rapport A 
l'axe Gz (la dimension de I z est celle d'une longueur A la 
puissance quatre) . 


2.1.2 - HypothAses 

En sus des hypotheses classiques : 

- les deformations restent petites. 


(H.l) 
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(H . 2) 


- chaque section droite reste plane au cours de la deforma- 
tion. 


on suppose que : 

(H.3) - l'axe Gz (oaralieie S l'axe de rotation K) est un axe 

principal d'inertie pour la section droite S. 


2.2 - CALCUL DU LAGRANGIEN DU SYSTEME 
2.2.1 - Preliminaires 

- Poutre non deformSe : le point G occupe la position Go 
de composantes (x,0,0) dans le repdre R. 

- Poutre deformSe : le point G occupe la position difinie 
par les composantes (x+u,v,0) dans le rep£re R. Si a(x) 
d^signe 1' angle que fait la tangente 3 la ligne moyenne de 
la poutre avec l'axe l, 1'hypothSse H.3 permet d'Scrire 


2.2.2 - Calcul de l'gnergie cinetique d'un element de la poutre 

Considerons une tranche eiementaire de poutre, comprise 
entre les sections droites voisines S(x) et S (x + dx) . 

L'energie cinetique de cette tranche s'ecrit comme somme : 

- de son energie cinetique de translation : 

dT t = I P' S dx 


- et de son energie cinetique de rotation qui (cf. hypothd 
se H.3) s' exprime sous la forme : 

dT r = 1 p T z w2 dx 

w designant la vitesse de rotation de cette tranche. 


2. 2. 2.1 - Ener 2 ie_cinetigue_de_ translation 

La vitesse absolue du point G est egale S la somme de sa 
vitesse d'entrainement V^ e par rapport au repSre inertiel I 
et de sa vitesse relative V Gr par rapport au repSre R. 



En d§signant par : 


to = 0 K 
e 


la vitesse de rotation d' entrainement du repdre R par rap- 
port au repSre I , il vient : 


V Ge = u. e A OG 

Projetons cette relation sur le repdre R. On obtient : 

-0 v 
0 (x + u) 
0 

' l 

De m§me, dans le rep§re R : 

u 
• 

v 
0 


0 


x + u 



* 



0 

OG 

V 

V Ge 

• 

0 


0 



'Gr 


II en risulte 1' expression de V G dans le repdre R 


-f 

V, 


u - Sv 
V + S (x + u) 
0 


et finalement : 


( 1 ) 



2. 2. 2. 2 - Ener 2 ie_cin§tigue_de_rotation 

Soit to r la 
repSre R : 


Soit to r la vitesse de rotation relative de S par rapport au 


w^. = a K 


I OS 



II vient : 


v-* 

03 


“e + “r 


soit : 


10 = 


(0 + a) K 


ce qui, puisque (cf. § 2.2.1) : 


a = v' 


fournit : 


9 • 


(0 = (9 + v') K 


Ainsi s 


( 2 ) 


dT r = 1 p I z (9 + v') 2 dx 


2.2.3 - Calcul de l'energie potentielle d'un Element de 
poutre 

Cette dnergie potentielle est egale a la somme de l'energie 
potentielle de traction-compression et de l'Snergie poten- 
tielle de flexion. 

2. 2. 3.1 - Energie_potentielle_de_ tract ion-compress ion 
Elle s'Scrit (cf. r6f§rence (1), tome 2, page 35) : 


1 N 2 

dV c = 2 ES dx 


avec : 

N : contrainte normale 

E : module d' Young du matSriau 

Comme (cf . r§f§rence (1) , tome 3, page 125) : 


N 
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c vs x wi 


il vient : 


(3) 


dV c = ES (u ' ) 1 dx 


2.2. 3.2 - Ene r 2 i§_go ten tielle_de_ flexion 

Elle s'^crit (cf. r§f§rence (1), tome 2, page 61) 

i M 2 
dV £ = ^ Z 


f 2 El 


oil est la contrainte de flexion, 
z 

Comme (cf . r§f§rence (1) , tome 3, page 125) : 


«z - «, It 


il vient : 


dV £ = 


1 

2 z 


El. (a ') 2 



ou encore (cf. § 2.2.1) 

(4) 

2.2.4 - Lagrangien du systSme 
Ce Lagrangien s'§crit sous la forme 

L = /t L dx 


avec : 


t = f - V 


et (cf. § 2.2.2) : 


t = ^ ps{[u- ev] ! + [v + e(x + u)] 2 } tipi (6 + v') 
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et (cf. § 2.2.3) : 


V = | ES (u') 2 + i EI z (v") 2 
D'oCl 1' expression de L : 


(5) 


L = i p S { [u - 9v] 2 + [v + 6 (x + u)] 2 } + i p I z (0 + V') 2 

- \ ES (u ') 2 - EI Z (v") 2 


2.3 - OBTENTION DES EQUATIONS DU MOUVEMENT 
2.3.1 - Formulation g§n 6 rale 
La poutre est soumise H un couple : 


M (t) = M (t) K 


au point 0 .‘ 

Les forces de pesanteur ne travaillant pas (car perpendicu- 
laires a tout dSplacement virtuel ( 6 u, 6 v, 50)), n'intervien- 
nent pas. 

Les equations du mouvement s'§crivent done sous la forme 
(cf . reference (2) , page 233) : 


(e.l) 


3L 3 ,3L . 

3u " 3x W 


A - 


3 x 2 9 u" 


3 ,3L. 

9t 3u 


3x3t 


A =o 
3u' 


(e. 2 ) 


3L 3 , 3L . 3 2 . 3Lv 

15v “ 3x '3v' ; 3x 2 '3v" 




3 , 3L , 

. + TTv'SV ' ~ ° 

3v 9x9t 3v' 


/_ o\ d r3L-i 3L _ 

<e - 3 > at [ ^J'3e- M 

avec les conditions aux limites : 

/ i i \ r3L 3 / 3L \ 3 , 3L x -i r 

(cl.l) tgp - ( 3 ^ 77 ) - ji (— )] «u = 0 pour x = 0 et x = L 

o \1 

/OB 



(cl. 2) 

3u" Su ' ’ 0 


pour x = 

(cl. 3) 

i-9L _ 3 / a l . a / a£, x -i r _ 

*-3v' 3x 3v" at ^ “ 

3v' 

0 

pour x = 

(cl. 4) 

- o 


pour x = 


2.3.2 - Formulation explicite 

Tous calculs faits, les relation 
sous la forme : 

(e.l) 

et (e. 2) 


(e' .1) pS 0 [v + §(x + u)] +A (ESu’) + p S (-u + 6v + §v) = C 

(e‘ .2) -pSMu-ev) - —(El v") - p S [v + e(x + u) + Su 

3x 2 2 


.± (P 
ax 

v 

Sachant que : 



la relation (e.3) s'§crit sous forme int§grale : 

(e ' • 3) [%] t = / t fc M(y) dy 

90 to 

avec : 

^7 = [j() L { I z + S [(x + u) 2 + v 2 ]} pdxl 6 

+ j[) L { I z V' + S [v (x + u) 


0 et x = L 
0 et x = L 

0 et x = L 
s 'expriment 

1 

I z v‘) =0 


- uv] } pdx 
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En ce qui concerne les conditions aux limites (cl.l) et 
(cl. 2), on obtient respectivement : 


3L 3 , 3L . 3 . 3L . 

“ 3x '3u"' " 3 1 


ESu' 


3L 

3u« 


0 


En ce qui concerne les conditions aux limites (cl. 3) et 
(cl. 4), on obtient respectivement : 


3L 3 / 3L i 

3v< ” 3x '3v«' 



(Jkj 

3v' 


= 4 (EI Z V“) 


- p I z (6 + v') 



-El, v" 
z 


a) conditions aux limites en x = 0 : 

En x = 0, la poutre est encastr§e. Elle satisfait done les 
conditions geom§triques : 

u(0,t) = 0 

v (0, t) = 0 

v' (0,t) = 0 

On a alors : 


6u | =0 6v | =0 <5v ' . =0 

I x=0 I x=0 I x=0 

et les conditions aux limites (cl.l) & (cl. 4) sont satis- 
faites. 


b) conditions aux limites en x = L : 

En x = L, la poutre est libre. Elle ne satisfait done aucune 
condition g£om§trique, ce qui permet de choisir des d^place- 
ments virtuels non nuls : 


6u | ^ 0 6v * ^ 0 

|x=L 


Sv' 


* 0 


no 


X=L 


X=L 


^z> 


Les conditions aux limites (cl.l) a (cl. 4) fournissent alors 
respectivement : 


ESu'i = 0 
I x=L 


u' (L, t) = 0 



v") 


- P I, 


(§ +?, )l|x= L 


0 


-El v" 
z 


x=L 


0 


V" (L,t) = 0 


III 


CONCLUSION 


Les Equations du mouvement d'une poutre flexible encastr€e 
a sa base sur un axe vertical en rotation s'gcrivent : 

. Squation de traction-compression 
(ESu 1 ) - pS\i}+pS [0 v + 2 6 v + e 2 (x + u)] = 0 
. Equation de flexion 

(p Z z v '> - rr < E v ”> - p s vj 

oX 

- p S [e (x + u) + 2 0U - fi 2 v] = 0 
. Equation de rotation d 1 ensemble 


3L 

ae t ‘ 

3L 

a Ht=t. 

II 

avec : 



& - i# 

30 L u 

+ s 

[(X + l 



+ 1 


Les conditions aux limites s'expriment sous la forme : 

. en x = 0 

u (0, t) = 0 
v (0, t) = 0 
v' (0,t) = 0 

• en x = L . 
u' (L,t) = 0 

(E I 2 v") - p I z (8 + V') =0 pour x = L 

1)2 


v" (L » t) = 0 


Flexible Beam 
Simulation 

by 

Shalom Fisher 
Tom Posbergh 
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Washington, D. C. 
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Flexible Beam Model. Nonlinear Slew Simulation 

• DISCOS simulation of nonlinear slew of SCOLE 
configuration with a flexible beam. 
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DISCOS CHARACTERISTICS 
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well as linear oscillations of 
system. 

External torques and forces can be 
added by the user. 



SCOUr CONFIGURATION 
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COMPARISON OF MObEL FlbELITY 
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Regulator Design 
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• The shuttle~ref lector-beam system is modelled by 
deterministic methods- 


Regulator for flexible 
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Closed loop damping and frequency shift 
Eigenvalues of closed loop system 

Defined modal frequency Closed loop eigenvalues 
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ISSUES IN MODELING THE SCOLE CONFIGURATION 

REVIEW HU DEVELOPMENT OF SYSTEM SOFTWARE FOR LSST 
DYNAMICS ANALYSIS. 

OPEN-LOOP SYSTEM DYNAMICS WITH RELATED STRUCTURAL 
ANALYSIS REPRESENTS FUNDAMENTAL STEP. 

FORMULATION METHODS. 

WHAT CAN WE LEARN ABOUT THE OPEN-LOOP SYSTEM? 



LSSI contract 

/operational 
IP- in progress 

Fig. 1 Development of system software for LSST 
dynamics analysis 
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FORMULATION TECHNIQUES 


(1) Eulerlan vs. La Grangian 

^ (2) Modeling of the Flexible Appendage (Mast) 
with Offset Inertial Masses at Both Ends 

(2a) Initial 2-D Vibration Analysis 
(1) P.D.E. approach 

(li) Finite element methods 

(ill) Treatment of boundary conditions 
Mast as cantilever 

Mast as a uniform beam with end bodies 
having inertia 

(iv) Separate treatment of lateral and torsional 
modes . 

(2b) 3-D Vibration Analysis 

D. K. Robertson, Sept. 1985 
Harris Corp. Dec. 198A 

(3) Linearization of 2-D Equations 

(A) Formulation of 3-D Equations Based on 3-D Vibration 
Analysis 
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II. A Development of the Two Dimensional Model - (Eulerian Moment 
Equations) 


The SCOLE system is assumed to be comprised of three main parts 
(Fig. 2.1): 

i) the Space Shuttle Orbiter with its center of mass located at 
point 0^; 

ii) the mast, treated as a 130 ft long beam, connected to the 
Shuttle at 0£ and to the reflector at 0^; 

iii) the reflector, considered to be a flat plate with its center 
of mass at 0^. 

The preliminary analysis presented here started before it was speci- 

Q 

fied that the interface point between the mast and the Shuttle is at 
0^. Therefore, in what follows, a position vector R^ appears which de- 
fines * where O 2 is the assumed interface point. 

In the following analysis, the angular momentum of the entire system 
is evaluated at point 0^ and the dynamics include the lateral displacements 
of the beam . 

II. A. 1 Angular Momentum of the Shuttle with Respect to Point 



O 
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Consider a point, P, of mass, dm, at an arbitrary position in the 
Shuttle such that O^P * r. The elemental angular momentum of the mass, 
dm, is given by: 

= Tx | 5? U. d " . = !*i cto k dm 

a rx^R-fc + Rtfcl/ * (to,- Jam (4.1) 

The total angular momentum for the Shuttle is obtained by integrating 
Eq. (2.1) over the entire mass of Shuttle as: 

f^ 0 = - Rt* f r*dro - RcOc't* C^dm + 

The first and second integrals appearing in the right side of Eq. (2.2) 
vanish because the center of mass of the Shuttle is at point 0^. 

Since r • j * 0, Eq. (2.2) takes the form: 

H */o, = («c-8)j - It 3^ CZ-3) 

where II- is the Inertia tensor of the Shuttle at point 0, and u)_ , » 

. ^ 1 R l /R 0 

(“ c - 9) j. 

II. A. 2 Angular Momentum of the Mast with Respect to Point 0^ 



1 36 





Consider here an element of the mast located at point, P^, with 
mass, dm. The elemental angular momentum of such an element Is given 
by= s jc X d dm 


(2.4) 


' — » « 

if one notes that 


(2.5) 


0P t r & + 


then, Eq. (2.4) may be expanded according to: 

/(' 

'V-m 


Hh/o, z * J ( r ° + JyLrdsfijtikk+jiifa ( 2 . 6 ) 


<M / 


f fU 


is expressed using the relationship between 


(2.7) 


the rate of change of a vector, w, in an inertial (R q ) and rotating 
(R^) frames, i.e. 

After substitution of Eq. (2.7) into Eq. (2.6) and integration term by 
term, one can develop: 

Hn/Of r A £ + /*»(**&) J °dr 

— J— <x*(u)b+<p)(£/wf&-6cnfl& + C +2> +8 -^y (2 ‘ 8) 

fit 

+ ^ + £ ' JfrnfiC - C cod -A + cjj + J 


Tf Qfat) sV> xO^uyuad - /f cc^T-f B/>»t/£z. + C ccd/QC + 


13 ? 


II. A. 3 Angular Momentum of the Rigid Reflector with Respect to 



Ox 


Let 0^ be the center of mass of the reflector, and 0^ the inter- 
face point between the reflector and the mast. The distance, X, between 
0^ and 0^ is constant since the reflector is assumed to be rigid, at least 
for this analysis. 

Let us now consider an element of mass, dm, of the reflector lo- 
cated at an arbitrary point, P 2 * The elemental angular momentum of that 
element of mass can be expressed as: 

dH. a 6 Tp,x i C0PJL dm 

r /o. dt UU 

0^P 2 and 0^2 can be expressed as: 
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Eq. (2.9) may be expanded according to 

d 


dH r/0/ = (W, *<k+x)U)*± (4+ *»*&**&$ fa 

fc+’&J /fa 


( 2 . 11 ) 


Once more, 
expressed using Eq. (2.7): 


is 




After substitution of Eq. (2.7) into Eq. (2.11) and integration term by 
term over the entire mass of the reflector, one arrives at 


Hr/ 0f - J Mr RR 2 oO c ' Cat * (*+&) f MrX 'RtA. * 


( 2 . 12 ) 

a 

where is the moment of inertia of the reflector about the j axis 


taken at point 0^. 



II. B.l Moment Equation 


The angular momentum of the entire system about 0^ is obtained 
by summing the angular momentum of each part about 0^, i.e. 


The moment equation 


M*/0 y r ^ u /o, 


(2.13) 



(2.14) 


where N is the sum of all the external torques, acting on the entire 
system, about an axis through point 0^. 

At this stage of the analysis, it is assumed that the center of 
mass of the Shuttle moves in a circular orbit, i.e. 


. *** A w* 

£L %/ — ~ & 

cbk JZ K (2.15) 


Taking into consideration the coincidence between points 0^ and 0^* 
Eq. (2.14) is expanded using once more Eq. (2.7) and the following 
result is obtained: 

&%h-t= ** - 

- §(Zu + + Mr(JC*t£)) 

- ( Z.r 4 Mr(X l +Zi)) 4 -L Mm U3<.R, 

4 (i+«)j Mm Ju> e .zj[&rt)£ 4 - Ifnliat+ty* 


14-2 



VJj 

_ Mr (ftR* cdt (Q+oS) + 3CR* °^)J^ 

^ Mm cDc. f (^ ^ ~ ^A/dm/SC 


— & O&jid -h C y&piA /%£ +£ c**& -+ 3- 

h M«* Cd Z yOr> Cwt+<t> )jl (A/bij$C+ C/faf)ftC - Bajygc 
4 2 ^ J-~ ^ A Ccnfti •/• B /fr*feC _ C anrfifoC 

— — A -f £ J 


Rsu£ 


OS s J- 6n (u>t+(p) (pC^) 
£ 

V = (_U>t+(p) (^({, J 


( 2 . 16 ) 
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FORMULATION TECHNIQUES 


(1) Euler ian vs. La Grangian 

(2) Modeling of the Flexible Appendage (Mast) 
with Offset Inertial Masses at Both Ends 

^ (2a) Initial 2-D Vibration Analysis 
(1) P.D.E. approach 

(11) Finite element methods 

(111) Treatment of boundary conditions 
Mast as cantilever 

Mast as a uniform beam with end bodies 
having inertia 

(iv) Separate treatment of lateral and torsional 
modes. 

(2b) 3-D Vibration Analysis 

D. K, Robertson, Sept. 1985 
Harris Corp. Dec. 1984 

(3) Linearization of 2-D Equations 

(4) Formulation of 3-D Equations Based on 3-D Vibration 
Analysis 



II. B.2 Expression for q 


In the moment equation, Eq. (2.16), one notices integrals involving 
q, the transverse displacement vector, and its first and and second deriva- 
tives with respect to time. It is, therefore, necessary to develop an ex- 
pression for q. 

II.B.2.i Relation between q(x,t) and y(x,t) 



Consider the beam in its deflected configuration, y (i ,t) is the 
deflection of the reflector-end of the mast at an arbitrary time^ t; 
y(x.,t), the deflection of an arbitrary point on the mast at the same 
time: . 

A 

From Fig. (2.1), * cos a (2.17) 

Assuming a small, tan a can be expressed as 



■x. 


(2.18) 


14-5 


(2.19) 


From Eq. (2.18) one derives 

C|(x,t) = x V.M - y(*>b) 

or 

<j(*> t) = <**- 


II. B. 2. 11 Evaluation of y(x,t) 

Assuming separability of the variables, the beam equation, 

+ 11 t ) _ 

JA" '»*•*' 


is solved to yield solutions of the form: 

3C*,fc) = ftt) 4>tx) 


where 

f(t) ■ E sinut + F cosut with u> ■ frequency of the vibration 
and <j»(x) » A cosSx + B singx + C coshBx + D sinhBx 

9 

When the following boundary conditions are assumed: 
a) y(0,t) ■ 0 ; b) y'(0,t) » 0 
c) El y'" (£,t) - -Mr yU.t); d) El y”U,t) = 0 


where Qna u 

< 1:1 

these can be expressed in the form: 

A + *58 so) 

Y A + fB so) 


U 

"it 

* 

<* 

1 






c 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 
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where M_ 

a ■ sinBi. - sinhBi - — r» B (cosBfc - coshBi) 

pA 

Mr 

6 « -cosBA - coshBi - — r> B (sinBfc - sinhBi) (2.26) 

p A 

Y ■ cosBA + coshBi 
a * sinBA + sinhBi 



Q 

For the SCOLE system, the following parameters have been supplied : 
pA 7 ■ 0.09556 slugs/ft 
El - 4. 0x10 7 lb-ft 2 
Mj. - (400/32.2) slugs 
l - 130 ft. 

For non-trivial solutions for A and B, det C must vanish. The 
values of jl for which det C ■ 0 are computed and substituted back into 
Eq. (2.28) to obtain the frequencies of the different vibrational modes 
(Table 2.1). 

The same values of B are substituted into $ (x) , (Eq. 2.23) J which is 
normalized with respect to its maximum value and the normalized mode 
shapes plotted (see Table 2.1 and Figs. 2.2 - 2.6). Note that the ranges 
of frequencies obtained in Table 2.1 a.re higher then those prevously pre- 
sented in the April 13, 1984 oral presentation due to previous inconsis- 
tencies in dimensional analysis of some physical units. 
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TABLE 2.1 


Values of 8 
for the First 

£ 

1.874599 

4.6929 

7.8519 

10.997 

14.1309 

17.276 

20.4229 

23.555 


and Natural Frequencies (HZ) 

8 In-Plane (Pitch) Bending Modes 

to (Hz) 

.677828 

4.245 

11.884 

23.3128 

38.4933 

57.5283 

80.4045 

106.958 
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a.cea 13.000 20.09a 39.009 £2.000 os. 000 73.000 91.000 10*. 000 117.030 130.000 

Fig. 2.4 



ISO 


Fig. 2.5 








II. C Frequencies of the Lateral Vibrational Modes when the SCOLE 

System is Modelled as a Free-Free Beam with End Bodies having 
Inertia 


The solution to the beam equation (2.21) is again considered 
and the following boundary conditions assumed: 

1. The shear force at either end is equal to the mass located at 
that end multiplied by the acceleration of the interface point at that 
end. 

This boundary condition combined with the equilibrium conditions 
yields 

_ £X Qjjffat ) _ Mi 

Oz 3 Qt*- 

at the Shuttle end, 

_ El - Ms - - Ms CL> Z yfajt) / 

at the reflector end 


PT yxfat)! _ Kir $&**}/ - _ Hr ti ) 1 U.so) 

ks/so * Iz-ljo 


where 


CD 3 , S 4 EJ. 

1 fA ’ 


1 52 



2. Next, expressing the equality between the moment at an end point 
and the product of the inertia of the mass at that end by the angular 


acceleration of the interface point results in: 

Ti BCzj t) = ET 

(vktAJL t) = ^ it (■** tl ^ 

<iz. 

At the Shuttle end: x * 0, this is expressed as: 



_ u> l Is 1 - EJ 7 l X(z,t )l (j.j/J 

Vz fxzO 9 Z 2 'Z=0 


the same boundary condition at the reflector end translates as: 

. co l If gflS&SV =- FI (J.u) 

^ x * Z -/£0 / z=JSO 


After performing the required differentiation of the assumed 
solution of the beam equation (Eq. 2.23), one arrives at the following 
system of four equations with 4 unknowns, A,B,C, and D, 


Eq. (2.29 A 
Jft' 

Eq. (2.30) =* 


+ 


8 + tiiA C - 2>-0 

JA> 
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Eq. (2.31) ap - A + Is 8 + C + 2sj£ J) s 0 (A-!S) 

■ I A' fA‘ 

and Eq. (2.32) 4> 

£ 2rJ$* Cfofot 1 A - [ lrj& * + /bu,$C 1 & 

JA* L fA> 

s. - If C + Tr /?> * Cad /%Cj - Q (£. 

$ A ^ A* 


Equations (2.33) - (2.36) can be recast in the matrix format as 



For non-trivial solution of <J>(x) (Eq. 2.23) the determinant of M(8) 
must be zero. A computer program was written, and the values of 8 , 
solutions of the nonlinear equation det [M(B)] - 0, obtained. 

These values of 8 were substituted into 



to derive the frequencies of the inplane and out -of -plane lateral 
vibrational modes. The results are given in Tables (2.2) and (2.3). 


IJ4- 



Table 2.2 Values of 3 and Natural 


Frequencies (H ) for the first 9 
z 

In-plane (Pitch) Bending Modes 


3 

“<V 

0.0097 

0.3065 

0.0310 

3.1308 

0.0549 

9.81922 

0.0789 

20.2809 

0.1030 

34.562 

0.1271 

52.6288 

0.1512 

74.4794 

0.1754 

100.229 

0.1995 

* 

129.664 


/ 



Table 2.3 Values of 8 and Natural 


Frequencies (H ) for the First 9 
z 

Out-Plane (Roll) Bending Modes 


8 

0.0103 

0.0310 

0.0549 

0.0789 

0.1030 

0.1271 

0.1512 

0.1754 

0.1995 


0.3456 

3.1308 

9.81922 

20.2809 

34.562 

52.6288 

74.4794 

100.229 

129.664 
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II. D. Derivation of the Frequencies of the Torsional Vibration, 

SCOLE Configuration. 

Assuming the mast to be a circular shaft, the torque at any point 
on the shaft is given by 

T- G1 2 b) 

^ oc 

where G is the modulus of rigidity and I the polar moment of inertia 

of the cross sectional area of the beam. This torque is opposed by the 
inertial torque 

J Qt Z 

where p is the density of the beam. For equilibrium, 

GJ ? &' *> _ I t Vj&tl -o fig.jg) 

lx J 7i z 


Assuming the separability of the variables. Equation (2.38) is solved to 
yield, solutions of the form 


where 


= fit) (p(x) 

f(t) - a cos (cut) + 8 sin (wt) 

§ (x) * A sin cu^p/g 1 x + B cos 


L H) 


J 
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Boundary Conditions 


Writing that the torque, T,at either and of the beam equals the 
moment of inertia times the angular acceleration of the interface 
point yields: 

GI 0?(Zytr J _ J; yjfjx.t) (S. *0) 
Oz Ot* 

Equation (2.40) along with the equilibrium of the shaft gives: 
for the Shuttle end: x * 0 

GlUtOZl/ =.- It / (£-■* 0 

QZ, IxzQ Izzq 

for the reflector end: x ■ l * 130 

_ GI ° / - - Ir ljS^Cx J b)l (J.4Z) 

Qz IZzliC tXztlC 

After substitution of equation (2.39) into equations (2.41) and 
(2.42), one arrives at: 

Eq. (2.41) ^ A GJ + 3 I S CO s O (4.43) 


Eq. (2.42) 



I SB 



Equations (2.43) and (2.44) can be recast in matrix format 



For non-trivial solution of equation (2.39) one must insure that the 
determinant [P(a>)] is equal to zero. 

The values of ui for which det [P (u>) 3 * 0 correspond to the 
frequencies of the torsional vibration. A computer program was written 
to solve this determinental equation and the frequencies for the tor- 
sional modes are listed in Table 2.4. 


1 5 ? 


Table 2.4 Values of Natural 
Frequencies (H z ) for the First 9 
Torsional Vibration Modes 


0.0305 

39.99 

79.98 

119.97 

157.97 
199.96 
239.55 
279.94 


319.939 



II. E Preliminary Calculation of the SCOLE Appendage Frequencies 
based on Finite Element Techniques 

For this application both the reflector and the mast are assumed 
to be a single flexible body. This body is considered to be comprised 
of two types of elements: (1) beam elements; and (2) triangular plate 

elements. The actual finite element model (FEM) is described as follows: 
Mass distribution 

Space Shuttle 6,366.46 slug 

Mast 12.42 slug 

Reflector 12.42 slug 

The masses of the reflector and the mast being so small (.39%) as compared 
with the mass of the orbiter, which in this analysis is ass um ed rigid, 
the system could be modelled as a cantilever beam (mast) with a mass with. 

(reflector) at its end. Also, the reflector in this section is 
going to be assumed flat with a constant thickness small as compared with 
its characteristic dimensions. 

The dynamics analysis of the STRUDL software package, which uses a 
physical analysis to solve the equations of dynamic equilibrium, is used 
to generate the eigenvalues, the frequencies, and the periods of the system. 


System Geometry (Model) 



The beam (mast) will be divided into 3 beam elements (each of 
43.33ft length) having a mass of 4.14 slug to be lumped at the 
ends of the elements. 








Results - Conclusions 


The following results have been obtained (Table 2.5). They show 

that the system is less stiff in this model as compared with previously 
8 

developed NASA and Howard University continuum models and also that re- 

9 

cently described by the Harris Corporation. 


Table 2.5 - Modal Frequencies (H ) 

Obtained by Implementing a FEM 
of the Preliminary Model of SCOLE 
(Poisson's ratio * 0.3 assumed) 


0.157 

0.275 

0.782 

1.083 

1.232 

1.386 

80.09 

107.24 

107.24 

265.99 

421.50 
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FORMULATION TECHNIQUES 


(1) Eulerian vs, La Grangian 

(2) Modeling of the Flexible Appendage (Mast) 
with Offset Inertial Masses at Both Ends 

(2a) Initial 2-D Vibration Analysis 
(1) P.D.E. approach 

(11) Finite element methods 

(ill) Treatment of boundary conditions 
Mast as cantilever 

Mast as a uniform beam with end bodies 
having Inertia 

(iv) Separate treatment of lateral and torsional 
modes . 

(2b) 3-D Vibration Analysis 

D.K. Robertson, Sept. 1985 
Harris Corp. Dec. 1984 

/ (3) Linearization of 2-D Equations 

(4) Formulation of 3-D Equations Based on 3-D Vibration 
Analysis 


t » 
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II. F 


Linear izaCion of the Equation of Motion— Floquet Analysis 

Let t the dimensionaless time be equal to a) t; — = u> — • 

c dt c dT ’ 



Dividing each term of Eq. (2.16) by M i yields 

- ® t + M * Lr / M ,^ i- Mr 3fx+ Mr } 

Mm 6 Mm J 

+|& ' *!) 101 W ^ +<*,)) 

~ Trlf ^ wltJn ^ t+ ^{( : ‘y Ml ,e» + Mr + Mr Vej 

+ ^ J - 03 Cuit+4) o»(wfc*<*,) 1 fie*. 

- dr %u5 t ] - ^ OM) «Kt) Mr X & <* * o ( 2 .«,) 

n* l* ) M«» t -t 

C ! {&*/M«|l4 ^ t + yU. +ylt j 

X 

where u= — ; X - j 
m 

c 2 

where ip (2) * A ^ jh* flJL + C 4 ^ 'SuUt 

- A * B 4 CyXuli -tD cw& ^ + B — ft 

C 3 ■ ( Xtr / M .e i +> ^l l )W + Hilti/Ji + ‘httjjfai 

where ^(i) ■ A* - & u»y$C 4 Cy>vulq ^ -t» 2> cpA 

^3 (i) - A C&* ft C B />K1 ^ - C OB5f^ ft/t - 1) / hu ll y^C — A 4 C 
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C 5 - yU,^R 
and Cg ■ 4*(.£) yU. ^ £l 
E q. (2.44) can be written as 

- C, 8 + O^Ca, J.(£ud^+<(>J) - 03 c C s 9 - + cO l Cj c®(o:fc-+ 4 >) 

« CJe.u ;c^ />^(wlr+$)<^(coM<t 9 - COCufe C £ - iiq* (jw-J^S’rO 

C2-47 

Introducing the dimensionaless time r - u» c t and dividing Eq. (2.47) 

by (o^ one arrives at 
7 c 

_ uCc, c, i (O’co.cujfc+^j) - u*?c s fl 4 - 3 ^ci„-i„)e- 

aZ dZ 

.f CO* CG^^tr+^O — COcUJ G4 Co (<ofc U)COcC^I^^O^-t <fj rO 

C2.4g) 

Introducing now the new parameter fl * — yields 

c 

_ G.W -t C* i i> <u»££LE+ 4 >)] _ C 5 sl B- 

at at 

4jd.*C J c*»(£it-Kj))- -0-t4 / 4 vm- Ta( rir -t- d 5 ^] 

A 

- £lC t J&* Inx +<$) - 3 (I„- I J3 ) & a O (2 AS) 


\G7 





FORMULATION TECHNIQUES 

(1) Eulerian vs. La Grangian 

(2) Modeling of the Flexible Appendage (Mast) 
with Offset Inertial Masses at Both Ends 

(2a) Initial 2-D Vibration Analysis 
(i) P.D.E. approach 

(li) Finite element methods 

(lii) Treatment of boundary conditions 
Mast as cantilever 

Mast as a uniform beam with end bodies 
having inertia 

(iv) Separate treatment of lateral and torsional 
modes. 

(2b) 3-D Vibration Analysis 

D.K. Robertson, Sept. 1985 
Harris Corp. Dec. 1984 

(3) Linearization of 2-D Equations 

^ (4) Formulation of 3-D Equations Based on 3-D Vibration 
Analysis 


i i 
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WHAT CAN WE LEARN ABOUT THE OPEN-LOOP SYSTEM? 

(1) Linearization of 2-D Equations (Torque-Free) 

(2) Stability Analysis - 

(2a) Assume appendage is vibrating at only one of 
its flexible modes 

(2b) System can be described by periodic coefficients 

(2c) Floquet analysis can be used to determine system 

parametric instabilities - auasi-analytic results 
obtained for cases of (i) No gravity-gradient, 
no offset of mast interface point on reflector, 

(ii) with offset of mast interface point, but 
not gravity-gradient. 

(2d) For general case, a numerical implementation of 
the Floquet analysis is required. 

(3) Relation between this system and other systems involving 

geometric offset 

(3a) The dynamics of orbiting tethered Platform systems 

2- D analysis, Stanley Woodard LSSI 

3- D analysis (in-progress). Fan Ruying - 

Visiting Scholar, Beijing Inst, of Control Engineering 

(3b) The dynamics of the Wrap-Rib Antenna system - any 
published results? 
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Parametric Study of the System 


Let us assume that the interface point between the reflector 
and the mast is at the center of mass of the reflector 
-*• X-0 -*■ X * 0 » C c - c, 

J 0 

Under this assumption , the equation becomes 

_ 6" + & d [&«» car-Kt;)! + n* £> 

H aC Gy 

— S± a /*** [ z int+<|>)] — i (, i»i -Th) 9 s o . 5 o) 

which yields the following first integral 

- 6’+ & [&M>c*it+(t>)J t /»«( at+4>) 

+ g c» fc u.ro 

This equation can be plotted in the phase plane (0^0) for different 
values of p and Q. 

Floquet Analysis 

The angular motion about an axis perpendicular to the orbit plane 
is described by: 

[- £.* * C» nx j Q 1 _ 6 

, (2.5-2; 
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This equation can be recast into the following matrix format 


9* 


ft 


r u 


P» 


( 2.5-3) 


where 


[PCr)J 


+ £> cciiir 
c, c, 


£ £§* c } •+ « (X»» ■“ Ijj 


1 


o 


Since P(r) is a matrix with periodic coefficients, the stability of the 
motion will be analyzed in what follows using the Floquet theorem. 


Case 1 . No gravity gradient. No offset 


^Cc) 

CZ(T)] 


£}> cro .Qt 


c I 


l 1 

- CP(t)] [Z(T)] 


- Q fa* rzx 

C| 


o 
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2. Floquet theorem 


Equation (2.5J) can be written as 

♦" + P(t) ♦ ■ 0 (4.3) 

where p(x) * - < is a periodic function * 

with period T * 2ir/n . 

By setting y x - $ ; y 2 = *' »y£ r one could 
recast equation (4.3) in the following state variable 
form. 



Assuming that the quantities Z 11 (T)r Z (T) r Z (t) r 

1 2 2 X 

and Z 22* t * are known ( wheE « th e Zij (T) are the 

elements of the matrix [Z (T)J which satisfies the 

matrix equation [Z (t)J '« [E(t) J [ Z (t) 3 and which, 

when t«0 , equals the identity matrixd ) Jone can 

derive the stability conditions applying the Floquet 
14 

theorem which states: 

(i) if for any j r 1 1 1 1 > l, (where || s^|| 
represents the modulus of the eigenvalue, s^ ), the 
zero solution of equation (4.3) is unstable; 

(ii) if for all j , | | s j | | < l, the zero 
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solution of equation (4.3) is asymptotically stable; 

(iii) if for all jrl|sj<! is not greater than 1, 
but for some, j, 1 1 s ^ 1 1« 1, and each Sj of modulus 1 is 
distinct from the others , then the zero solution is 
stable but not asymptotically stable. 

(iv) if for all j, || s jJ| is not greater than 

1 , and there exist some with 1 1 s ^ | J « 1 in 

multiplicity n , the zero solution of equation (4.3) 
is unstable unless n - k where the s^ are the 
eigenvalues of I Z (T)],s n » repeated eigenvalues, and 
k - nullity of [ Z (T)-s n U -(order of I Z(Tl-s n I] 
minus rank of [Z (T) -s I]>. 



z u 

m 

”u z u 

+ 

p 12*21 

(1) 


Z 12 

- 

hl Z 12 

f 

P 12 z 22 

(2) 


Z 21 

at 

P 21 Z 11 

+ 

P 22*21 

(3) which becomes 

Z 21 * Z 11 slnce »2X 






1 and p 22 - 

0 

^22 

m 

P 21 Z 12 

+ 

P 22*22 

(4) which becomes 

*22 “ *12 

from (3) Z 21 * Z 

11 

substituted 

into 1 yields 



*21 " p ll*21 + P 12 Z 21 

Similarily from (4) Z 22 - Z 12 which substituted into (2) yields 

*22 “ P 11 Z 22 + p 12*22 
If one notices that p^ 2 ■ ~ p^ 

then Z 21 - P U Z 21 + Pn Z 2l " ft (p ll*21 ) 

and Z 22 - p u z 22 + P U Z 22 - fj- (p u z 22 ) 

These two last equations are integrated and the following result 
for Z 2 ± and Z 22 obtained 

*21 “ P U Z 21 + *1 

*22 " p 11*22 + K 2 
but from (3)^ Z 21 » Z^Ct) and 

from (4)^ Z 22 - Z 12 (t) 

therefore, Z* 21 (0) - Z^CO) - 1 - P u (0)Z 21 (0) + K. 
or for <> - 0 1 ■ K 1 since Z 21 (0) - 0 

-* K x - 1 

*22 (0) * Z 12 (0) * 0 * p ll (0) Z 22 (0) + K 2 
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or for $ ■ 0 


— - - Kj since CO) « 1 


Z 21 “ P 11 Z 21 + 1 


Z 22 " P 11 Z 22 ~ -T- 


Solution of the linear first order equation 


^22 

d T " P U Z 22 
dZ, 


(1) 


The presence of 2 2 and in the equation suggests a pro- 

dr ^ 

duct of the type ♦ C‘ c ) Z 22^ 


dr 


dr 


( 2 ) 


Multiplying (1) by 4»(r) yields 


$ at* _ ~ "|f ^ 


(3) 


which can become 


A(fe) = - & <f 

if one can find $(r) (the integrating factor) such that 

sji =-4fu 

dr 


(4) 


(5) 
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Z = f- 


& CZD^lT etc 
C| 


JU<bl£) * - S3L + K &r 

c s a 

<$tz) * «j> |- £& Aw* AX Hr VC J 


from 


4- t<K* ) 


Si 4^ 

C| 


one arrives at 


4** - J"-§f 4 Ct>dl 

" = J- f- & d>fe) dt 

<PW J ^1 

«?«cc)r e*f> [rk^ nr * k ] j _ ff ( ex f[^ /a ^ ric+k ^ c j 

According to Taylor’s series development of a function 

* K «f [-a ^itt] •«"{•-&* 4 0 - 

which is integrated term by term to give 



since 


Lc)-l*> K, s - |i 




Sj- ewfit. 9v^ flr ? f. Si + c - Si t z + 

Cj ‘ ^C( J Ci, Ci i- 
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cos J2t 


Solution of “ ^11^21 + ^ w ^ ere ■ — 

=> = 4 ^ ~^z\ 

&Z Az Az 


(i) 


- 4>> (I .?*, *4’ 

dt 


( 1 ) - ( 2 ) ^ 


rf = ' 


( 2 ) 


==> JL 4>’cx) = - Jkk. /W nx tK } «r 

X2* Cf 

4’(C)r -ex/jo Pm. D.Z •+ k J j 

|L » 4> J *> 4>*&. -J4> J at 

°£ii = .L ^clx - fh*Clz- ^c^cll 


According to Taylor's series 


4’fc)= 


/0O 



which is integrated term by term to yield 


• «r {' -M c ‘- $ v ’ *■’ ] 

2. 2.j fe) r O ^ - O 

CO =■ 


(? ) -- **r [<% * ar] [ f * - §f/ - J( to ^ -> ) 


It can easily be verified that 
and finally 

ku.s ^Lcj 


z '^°> = e'flSi ^ 0 ! 5 i 


V 


(tonr -i _)+ (l|-cmnt ex.|>|^. 
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With the use of a computer program, the eigenvalues of the lZ(t) ] 
matrix are computed for t ■ a period and their modulus compared with 
1 to determine the values of the parameters for which the system is 
stable. The results of such parametric study are shown in the 
following stability diagram. Fig. 2.7. The large number of unstable 
points in the parametric space (Q, MU) are thought to be attributed 
to the absence of the gravity-gradient torque in the model. Future 
plans call for the extension of the Floquet analysis for the cases 
where a non-zero reflector attachment offset is considered and also 
where both a non-zero offset and the effects of gravity -grad lent are 
included . 
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WHAT CAN WE LEARN ABOUT THE OPEN-LOOP SYSTEM? 


(1) Linearization of 2-D Equations (Torque-Free) 

(2) Stability Analysis - 

(2a) Assume appendage is vibrating at only one of 
its flexible modes 

(2b) System can be described by periodic coefficients 

(2c) Floquet analysis can be used to determine system 

parametric instabilities - quasi -ana lytic results 
obtained for cases of (i) No gravity-gradient, 
no offset of mast interface point on reflector, 

(ii) with offset of mast interface point, but 
not gravity-gradient. 

(2d) For general case, a numerical implementation of 
the Floquet analysis is required. 

(3) Relation between this system and other systems involving 

geometric offset 

i/ (3a) The dynamics of orbiting tethered Platform systems 

2- D analysis, Stanley Woodard LSSI 

3- D analysis (in-progress). Fan Ruying - 

Visiting Scholar, Beijing Inst, of Control Engineering 

(3b) The dynamics of the Wrap-Rib Antenna system - any 
published results? 
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SEPARATED LINEARIZED EQUATION WHEN 



WHAT CAN WE LEARN ABOUT THE OPEN- LOOP SYSTEM? 


(1) Linearization of 2-D Equations (Torque-Free) 

(2) Stability Analysis - 

(2a) Assume appendage is vibrating at only one of 
its flexible modes 

(2b) System can be described by periodic coefficients 

(2c) Floquet analysis can be used to determine system 

parametric instabilities - quasi -analytic results 
obtained for cases of (i) No gravity-gradient, 
no offset of mast interface point on reflector, 

(li) with offset of mast interface point, but 
not gravity-gradient. 

(2d) For general case, a numerical Implementation of 
the Floquet analysis is required. 

(3) Relation between this system and other systems involving 

geometric offset 

(3a) The dynamics of orbiting tethered Platform systems 

2- D analysis, Stanley Woodard LSSI 

3- D analysis (in-progress). Fan Ruying - 

Visiting Scholar, Beijing Inst, of Control Engineering 

/ (3b) The dynamics of the Wrap-Rib Antenna system - any 
published results? 
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ISSUES IN CONTROLLING THE SCOLE CONFIGURATION 

SHOULD CONTROLS ANALYSIS PROCEED IN TWO STEPS, 
FIRST USING A 2-D MODEL? 

WHAT IS THE EFFECT OF A TIME DELAY IN THE CONTROL? 



Systems with delay in Control 


Stability: 


In Literature: 

X(t) = AX ( t ) + BX(t-x ) 


A is a stable matrix 
stability of combined system is 
analyzed as a function of stability 
parameters of A and B matrices 


Control : 


In Literature: 


X(t) = AX ( t ) + BX ( t-T ) + CU( t ) 
U(t) = -A X(t) - B X ( t-x ) 
or 

X(t) = AX ( t ) + BU ( t-T ) 

U(t) = KX( t+T ) 
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Present Problem: 


X(t) = AX( t ) + BU( t ) 
U(t) = KX ( t-T ) 


A is marginally stable 
Design K such that controlled 
plant is stable. 



Con fret lane SirucJurt* 

w c-Jfi o tobuj os* Cpv\/tr^C • 


Xk> A>8h- BUtt) 

u(t)= kxet-T) 


k* ^JUjdv ’1 Kj 6&‘ oJ>ovt S*jtfen\. 


~ it StoJffe- 
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ORIGINAL PAGE IS 
OC POOR QUALITY 


B. Effect of Delay on Control Syst en Stability 

the control loo of the form U » -U Is designed for e systan 

rqn^tW . 

X = A/.+ BU <»•« 


sod s delay x is incorporated into the state variable, X, and its 
effect on control systea stability is analysed onnerlcally. 

To start with, a seco n d order differential equation representing 


the vibration of a structure in a particular node is considered. 

The syst an equation ne ed for wsrlc si analysis is given as 

d X 1 

jTi -+- to X = H <»•*> 

The control 0 is a rate feedback given by 

a 

U- - 2. S w gr| <*-3) 

and the closed loop equation with control is written as: 


dx 

dt a - 


+ 2^1) 




The miner leal values selected for numerical simulation are 
arbitrary and are: 


m ■ 

6.0 

C - 

0.5 

x(o) - 

0.5 

x(o> - 

0.0 


(B.4) 
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Transiant BetpooM of Second Oedar Systea without rdae Delay. 
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Figure 5*2 Transient ieapmac Ox Second Cedes System witu Time vela/ 
t ■ 0.1 secs. 
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B.3 Transient Response of Second Order System with Time Delay 
T *■ 0.15 secs. 
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1. INTRODUCTION 


In recent years much attention in the controls community has centered 
on problems posed by large space structures. These problems include the 
need to provide active-damping structural control for many, often densely 
packed, modal frequencies, and the need for distribution of actuation and 
sensing equipment over and throughout the large, lightweight structure. In 
the analysis and control synthesis of such problems, integral parts are 
played by the level of accuracy of the system model and the robustness of 
the controller. Properly done, the control design should provide damping 
to the vibrational modes within the controller bandwidth. 

This report accounts the progress made on control experiments for 
NASA's flexible grid experimental apparatus [1,2]. The grid is a 7-ft by 
10-ft lattice constructed of overlayed aluminum bars of rectangular cross- 
section 2-in by 1/8-in, centered at one-foot intervals, and is suspended on 
cables at two locations across the top bar. This design admits appreciable 
low-frequency structural dynamics, and allows for implementation of distri- 
buted computing components, inertial sensors, and actuation devices. Figure 
1 depicts the apparatus. Instrumentation includes non-contacting displace- 
ment sensors mounted on a separate rigid structure behind the grid, six rate 
gyros, and six inertia wheel actuators allowing application of torques up to 
20 oz.-in. 

A finite-element analysis of the grid provides the model for control 
system design and simulation. The motions of the grid perpendicular to its 
plane are of interest in this study, so that in the modeling analysis 88 
nodes are utilized with four degrees of freedom. 

The control strategy for this study involves a decentralized model 
reference adaptive approach using a variable structure control [3]. Local 
models are formulated based on desired damping and response time in a 
model-following scheme for various modal configurations. Variable structure 
controllers are then designed employing co-located angular rate and position 
feedback. In this scheme local control forces the system to move on a local 
sliding mode in some local error space. An important feature of this 
approach is that the local subsystem is made insensitive to dynamical 
interactions with other subsystems once the sliding surface is reached. 
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2. PROBLEM STATEMENT AND PROGRESS REPORT 


The overall experimental process for purposes of this study consists of 
three basic stages. First, the finite-element analysis is done for the 
grid apparatus on the Cyber 175 network computing system at NASA's Langley 
Research Center. The vehicle by which this is done is the SPAR analysis 

software package. Next, data generated by SPAR is used in a modal dynamics 
simulator, either on the Cyber 175 system or, as in the case of the results 
reported on here, on the Electrical Engineering Department's VAX 11/785 
system at The Ohio State University. It is within this simulation software 
that the control algorithm is implemented. The final stage of the process 
involves actual on-site testing of the design. The overall process is 
depicted in Figure 2. 

Modeling 

Assume that a general mathematical description of the grid structure 
takes the form 

MX + KX = F , 

where the matrix X contains incremental displacement variables for the 96 
grid points (88 for the grid, 8 for the cables) in each of the 6 degrees of 
freedom (3 axes, rotation about each axis). As noted above, typical 
modeling exercises involve only four degrees of freedom, where rotation 
about the axis perpendicular to the plane and translation along the horizon- 
tal grid axis are constrained. Also, M and K are the mass and stiffness 
matrices of the structure, respectively, and F is the vector of forces 
used in control. Since applications and environments of large space 
structures dictate the lack of any appreciable natural damping, the model 
excludes any damping. By employing the unitary transformation 

X = W , 

where ♦ is the mode shape matrix and W contains generalized coordin- 
ate modal displacements, a set of uncoupled equations results, namely, 

(# t M*)V + ($ t K*)V = * t F = u . 

In this expression, $ T M$ is the diagonal modal mass matrix, ♦*!($ is the 
diagonal stiffness matrix, and u represents the generalized forces 
(control inputs). In a truncation approximation from the finite-element 
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procedure, n modes are retained, implying that * T M* and * T K* are 
nxn in size. In a decentralized setting, 



where w i represents the i-th modal frequency. The development of the 
control algorithm relative to this viewpoint is given in the next section. 

The modal coefficients for this analysis (at gyro and actuator loca- 
tions) are taken from the SPAR program output, for a model retaining the 
first 10 modes. 

Software and Simulation Testing 

Simulation tests reported on in this study are all products of software 
implemented on the VAX 11/785 in the EE Dept, at OSU. The first stage of 
the code serves the purpose of transfering the modal coefficients from the 
SPAR output data set and formulating the proper mode shape matrix. Initial 
tests use one-mode models, resulting in a two-state grid model, and later 
simulations use the full ten-mode model, resulting in a twenty-state 
system. Next in the simulation code is the measurement, or sensor subrou- 
tine, from which the rate gyro feedback is obtained. Angular displacement 
is computed from these measurements with a simple numerical integration 
routine. After initial conditions are applied, integration is completed 
prior to updating of the computed control law. 

The grid is initially perturbed in the software by way of initial 
conditions on the angular displacements and their velocities; future 
studies will involve a sinusoidal excitation of the grid before control is 
applied. These initial conditions are introduced through the local refer- 
ence models, which are in all cases simple two-state (angle and angular 
velocity) systems. Various other segments were added to the code as needed, 
such as feedback of Kaman probe output and clipping constraints on the 
magnitudes of the applied torques. Primary output of the software is in the 
form of time-scale plots of linear displacement, rate gyro output, estimates 
of the angular displacements, and the applied torques, all at various grid 
node-points. Discussion of several tests is included in the next section. 
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3. DETAILS OF RECENT RESULTS 


Algorithm Development 

The Decentralized Model Reference Adaptive Controller with Variable 
Structure Control has evolved through four stages, the last stage specifi- 
cally for this application; 

(1) Variable Structure Control (Utkin, Itkis and others in USSR) 

Basically, the standard variable structure controller first drives a 
system's state trajectory to a given plane in state space and then, ideally, 
slides along this plane to the origin. The system can be shown to have 
excellent sensitivity properties while in the sliding mode. The control 
algorithm essentially checks on which side of the plane the states are and 
varies the feedback structure accordingly to orient the trajectories towards 
the plane. In actual implementation, the trajectories may "chatter" while 
sliding to the origin along the plane. 

(2) Model Reference Adaptive Control Using Variable Structure 
Feedback (Young) 

D. Young was one of the first to use varable structure controllers in 
driving the error signal to zero while doing model reference adaptive 
control. The approach retains all the advantages and disadvantages of 
variable structure controllers. 

(3) Variable Structure Controllers for Decentralized Model Reference 
Adaptive Control for Interconnected Systems (Ozguner, Morgan, 
Al-Abbass) 

A recent development has been the application of variable structure 
controllers for model reference adaptive control of interconnected systems 
with local state information availability. Ve summarize this approach 
briefly in the following. 

Consider the system 

N 

S i s *i = Vi + B i u i + EAijXj 

j*l 
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y* = d.x, , 


for i s 1, 2, N , where 


- „ni 

x. e R , 


A. e , 


y. e R J 


B. E R n i xl 


n i £ R 


D. £ R lxn i 


Given the model, the problem is to design a decentralized adaptive 
controller such that the states of each local subsystem are regulated to 
zero or track the states of a reference model. Each local controller is 
dependent only on the local subsystem, and is not allowed to communicate 
with the other local controllers. The only information provided for the 
local controller is the upper bound on the size of the dynamic interactions 
[ EA. .x . ] 

Let the local reference model for the i-th subsystem be given as 

* A 

x. = Ax. + B.r. 

1 11 11 

y. = x. , 

with x i in R n i , r j . the reference input, B A in R n i x * , and A i in 

R n i M i . Furthermore, let 

* t = C.e. , 

1 xn 

for C. in R i , specify the sliding surface, where 


e i - x i " x i 


and 


e. = A.e. + (A. - A.)x. + B.r. - B.u, - • 

3-1 

jHi 

The control lav is obtained from this last expression as 


•u, - *X * . Kjr, * , 


224 


where K r in r\ K p in R^* 0 *, and K # in R^^i can be specified [3] in 
regions of the state space, and where & i is a constant picked according 
to the norm of the interactions. The elements R r , K , R # , and S. are a 
function of the sliding surface and the coefficients of system and refer- 
ence model state equations. 

Another development used in the above is the idea of time varying 
sliding surfaces, which aide in smoothing of the chattering. 

(4) Multi-modeling for Decentralized Model Reference Adaptive Control 
using Variable Structure Controllers. 

Two developments were required to use the above approach for the 
control of flexible structures. The first was the incorporation of output 
instead of state feedback into the interconnected system model. The second 
and more crucial development was the expansion of the essentially general 
decentralized structure with multimodeling into an interconnected system. 

As an illustration, consider the two channel case 

x = Ax + BjUj + B 2 u 2 
Vi = D,x 
y 2 = D 2 x . 

We would like to consider two separate systems, 

{ A , B x , D x } , { A , B 2 , D 2 } 

i 

for designing the controllers. The system model used is 


r a * i 


r B x °1 


[ D ! °1 

L * A J 

9 

L o b 2 J 

9 

L o d 2 J 


where the "interconnection matrices" are not specified. This type of 
expansion has been previously utilized by Siljak in a different context. In 
the present context the specification of the fictitious interconnections are 
not required since the only information needed is an extra dominant term in 
the local control to suppress interactions. 
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Simulation Test Results 


In the following pages the results for two of several simulations 
testing the decentralized control algorithm are presented. For each case 
the group of output plots is preceded by a short summary depicting the 
number and location of inputs and outputs (for the numbering of the nodes, 
see Figure 1), program parameters (such as initial conditions, algorithm 
scale factors, and so on), and the structure of the reference model. Both 
simulations shown here are for the full 10-mode grid model. Note that the 
second of these tests represents results of sampled systems. 
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SIMULATION TEST #4 

10 modes 

6 inputs, 6 outputs, 

* co-located rate gyros and actuators 
located at nodes 25, 31, 52, 61 

* Kaman probe 2 with actuator 1 (output feedback) 
Kaman probe 7 with actuator 4 (output feedback) 

Program parameters: (w (0) and w.(0) not listed are 0) 

w 1 (0) = 0.1 1 

w 3 (0) = 0.1 
w 5 (0) = 0.1 
w 7 (0) = 0.1 
v 9 (0) = 0,1 
w 2 (0) = 0.2 
w 4 (0) = 0.2 
w 6 (0) = 0.2 
w (0) = 0.2 
MAX = 1 


Reference model eigenvalues: 

= -0.5 ; 


Full model (ten modes) 
output feedback, and reference model has 
critically damped modes. 

*** NOTE *** 

The applied torques are clipped (constrained) 
as 


X 2 = -0.5 

is used with position 


|applied torques | ^ 1.25 in-# 
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SIMULATION TEST #6 
10 modes 


6 inputs, 6 outputs, 

* co-located rate gyros and actuators 
located at nodes 25, 31, 52, 61 

* Kaman probe 2 with actuator 1 (output feedback) 
Kaman probe 7 with actuator 4 (output feedback) 

Program parameters: (w (0) and w.(0) not listed are 0) 

w 1 (0) = 0.1 1 

w 3 (0) = 0.1 
w 5 (0) = 0.1 
w 7 (0) = 0.1 
v 9 (0) = 0.1 
w 2 (0) = 0.2 
w 4 (0) = 0.2 
w g (0) = 0.2 
w g (0) = 0.2 
MAX = 1 

<V a 5 = 2 

Reference model eigenvalues: 

\ = -0.5 ; X 2 = -0.5 

Comments: 

Again, full model (ten modes) is used with 
position output feedback. 

*** NOTE *** 

This test represents the sampled 
version of the algorithm; sampling 
period here is 32 ms. 

Again, the applied torques are clipped 
(constrained) as 

|applied torques | « 1.25 in-# 
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IV. REMAINING WORK AND SUGGESTIONS FOR FUTURE STUDIES 


At the time of writing the present report a number of tests were 
completed on the ten-mode simulation model of the grid. In considering 
these results one can express confidence in the successful outcome of the 
application onto the real system. Both on-site and off-site (remote) 
experiments are to be performed during which the feedback control algorithms 
will be implemented on the Charles River Microcomputer system connected to 
the sensing and actuation on the grid. The gain values which gave satisfac- 
tory results on the simulation model are to be tried. 

A number of follow up studies can be envisaged using the software and 
hardware that is now available, and particularly using the expertise 
obtained. These include both theoretical studies which are required to 
answer some natural questions that arose during the tenure of the present 
project, and practical studies to test out some other reasonable algorithms 
on the grid problem and the results of the above mentioned theoretical 
studies. 

Specifically, further work is required to understand the implications 
of the multimodeling approach used in the present application. This model 
is especially suitable for local reduced order modeling. That is, one could 
try using only some of the modes to be controlled in one channel's model and 
some in the other, while some have to appear in both. The stabilizability 
implications and the spillover effects in using such an approach must be 
analyzed. An algorithm should be obtained for generating such different 
multiple models. 
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SUMMARY 

The problem of controlling large, flexible space systems has been the 
subject of considerable research. Many approaches to control system 
synthesis have been evaluated using computer simulation. In several cases, 
ground experiments have also been used to validate system performance under 
more realistic conditions. There remains a need, however, to test 
additional control laws for flexible spacecraft and to directly compare 
competing design techniques. In this paper an NASA program is discussed 
which has been initiated to make direct comparisons of control laws for, 
first, a mathematical problem, then an experimental test article is being 
assembled under the cognizance of the Spacecraft Control Branch at the NASA 
Langley Research Center with the advice and counsel of the IEEE Subcom- 
mittee on Large Space Structures. The physical apparatus will consist of a 
softly supported dynamic model of an antenna attached to the Shuttle by a 
flexible beam. The control objective will include the task of directing 
the line-of-sight of the Shuttle/antenna configuration toward a fixed 
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target, under conditions of noisy data, limited control authority and 
random disturbances. The open competition started in the early part 

of 1984. Interested researchers are provided information Intended to 
facilitate the analysis and control synthesis tasks. A workshop is planned 
for early December at the NASA Langley Research Center to discuss and 
compare results. 


INTRODUCTION 

Many future spacecraft will be large and consequently quite flexible. 
As the size of antennae is increased, the frequencies of the first flex- 
ible modes will decrease and overlap the pointing system bandwidth. It 

will no longer be possible to use low gain systems with simple notch 
filters to provide the required control performance. Multiple sensors and 
actuators, and sophisticated control laws will be necessary to ensure 
stability, reliability and the pointing accuracy required for large, 
flexible spacecraft. 

Control of such spacecraft has been studied with regard given to 

modeling, order reduction, fault management, stability and dynamic system 
performance. Numerous example applications have been used to demonstrate 
specific approaches to pertinent control problems. Both computer simula- 
tions and laboratory experiment results have been offered as evidence of 
the validity of the approaches to control large, flexible spacecraft. 

Concerns remain, however, because of the chronic difficulties in control- 
ling these lightly damped large-scale systems* Because of these concerns 
and because of the desire to offer a means of comparing technical 

approaches directly, an NASA/IEEE Design Challenge is being offered. An 


238 


experimental test article is being assembled under the cognizance of the 
Spacecraft Control Branch at the NASA Langley Research Center with the 
advice and counsel of the IEEE (COLSS) Subcommittee on Large Space 
Structures. This Spacecraft Control Laboratory Experiment (SCOLE) will 
serve as the focus of a design challenge for the purpose of comparing 
directly different approaches to control synthesis, modeling, order 
reduction, state estimation and system identification. 

The configuration of the SCOLE will represent a large antenna attached 
to the Space Shuttle orbiter by a flexible beam. This configuration was 
chosen because of its similarity to proposed space flight experiments and 
proposed space-based antenna systems. This paper will discuss the "Design 
Challenge" in terms of both a mathematical problem and a physical experi- 
mental apparatus. The SCOLE program is not part of any flight program. 

SYMBOLS 

a acceleration vector ft/sec^ 

A beam cross section area 

c observation matrix 

d noise contaminating direction cosine matrix measurements 

e line-of-sight error 

E modulus of elasticity 

f concentrated force expressions 

F/ 4 force vector 

g concentrated moment expressions 

GI torsional rigidity 

I moment of inertia matrix for entire Shuttle/antenna configuration 
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moment of inertia matrix, Shuttle body 

moment of inertia matrix, reflector body 

beam cross section moment of inertia, roll bending 

beam cross section moment of inertia, pitch bending 

beam polar moment of inertia, yaw torsion 

length of the reflector mast, beam 

control moment applied to the Shuttle body 

control moment applied to the reflector body 

disturbance moment applied to the Shuttle body 

mass of entire Shuttle/antenna configuration 

mass of Shuttle body 

mass of reflector body 

mass density of beam 

beam position variable 

direction cosine matrix, Shuttle body () . = T, ()„. , , . 

* J earth 1 Shuttle body 

direction cosine matrix, reflector body () « T,() __ 

’ 1 earth 4 reflector 

body 


inertial velocity. Shuttle body 

inertial velocity, reflector body 

lateral deflection of beam bending in y-z plane 

lateral deflection of beam bending in x-z plane 

angular deflection of beam twisting about z axis 

position variables 

displacement of proof-mass actuator 

line-of-sight pointing requirement 

noise contaminating angular velocity measurements 
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0,$,T pitch, roll, heading 
C damping ratio 

T noise contaminating acceleration measurements 

<0j angular velocity of Shuttle body 

<04 angular velocity of reflector body 

DISCUSSION 

The objective of the NASA-IEEE Design Challenge concerning the control 
of flexible spacecraft is to promote direct comparison of different 
approaches to control, state estimation and systems identification. The 
design challenge has principal parts, the first using a mathematical model, 
and the second using laboratory experimental apparatus. The specific parts 
of the Spacecraft Control Laboratory Experiment (SCOLE) program will be 
discussed in detail. 


Control Objectives 

The primary control task is to rapidly slew or change the line-of- 
sight of an antenna attached to the space Shuttle orbiter, and to settle or 
damp the structural vibrations to the degree required for precise pointing 
of the antenna. The objective will be to minimize the time required to 
slew and settle, until the antenna line-of-sight remains within the 
angle 6 . A secondary control task is to change direction during the 
’’on-target” phase to prepare for the next slew maneuver. The objective is 
to change attitude and stabilize as quickly as possible, while keeping the 
line-of-sight error less than 6. 
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Math Model Dynamics 


The initial phase of the design challenge will use a mathematical 
model of the Shuttle orblter/antenna configuration. It Is necessary to 
obtain a balance, of course, between complex formulations which might be 
more accurate and simplified formulations which ease the burden of 
analysis. 

The dynamics are described by a distributed parameter beam equation 
with rigid bodies, each having mass and inertia at either end. One body 
represents Space Shuttle orbiter; the other body is the antenna reflector. 
The equations for the structural dynamics and Shuttle motion are formed by 
adding to the rigid-body equations of motion, beam-bending and torsion 
equations. The boundary conditions at the ends of the beam contain the 
forces and moments of the rigid Shuttle and reflector bodies. The 
nonlinear kinetmatlcs couples the otherwise uncoupled beam equations. 
Additional terms represent the action of two, 2-axis proof-mass actuators 
at locations on the beam chosen by the designer. 

The rigid-body equations of motion for the Shuttle body are given by: 

- - I I 1 <Vi“x + M i + "d * Yi> 


S 1 
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Similarly, for the reflector body. 
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The direction cosine matrices defining the attitudes of the Shuttle and 
reflector bodies are given by: 
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The direction cosine matrices defining the attitudes of the Shuttle and the 
reflector bodies are related to the beam end conditions. 
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The equations of motion for the flexible beam-like truss connecting the 
reflector and Shuttle bodies consist of standard beam bending and torsion 
partial differential equations with energy dissapative terms which enable 
damped modes with constant characteristics for fixed, though dynamic, end 
conditions. The system of equations can be viewed as driven by changing 
end conditions and forces applied at the locations of the proof-mass 
actuators. 
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{PROOF-MASS ACTUATOR} 


{REFLECTOR BODY FORCE} 


{SHUTTLE BODY FORCE} 


{PROOF-MASS ACTUATOR FORCE} 


{PROOF-MASS ACTUATOR FORCE} 


{REFLECTOR BODY FORCE} 
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The angular velocity of the reflector body is related to the Shuttle body 
by: 
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The line-of-sight error described in figure 2 is affected by both the 
pointing error of the Shuttle body and the misalignment of the reflector 
due to the deflection of the beam supporting the reflector. The line-of- 
sight is defined by a ray from the feed which is reflected at the center of 
the reflector. Its direction in the Shuttle body coordinates is given by: 


- R R + R r t2 

[ r I <r r - v • v 

IK- “f- 2 

1 1 

• R a| 


where : 


Rp is the feed location (3.75, 0, 0) 

r R is the location of the center of the reflector (18.75, -32.5, 
-130) 

r A is a unit vector in the direction of the reflector axis in 
Shuttle body coordinates 


The vector R^ can be related to the direction cosine attitude matrices 
for the Shuttle body, Tj , and the reflector body, T4, by 



The relative alignment of the reflector to the Shuttle body is given by 
T 

TjT^ which is a function of the structural deformations of the beam. 
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The line-of-sight error, e, Is the angular difference between the 
target direction, given by the unit vector, D x , and the line-of-sight 
direction in Earth axes, T^los* 

« ■* ARCSIN |d t x Tj E los | or ARCSIN 

Computer programs are available which generate time histories of the 
rigid body and the mode shapes and frequencies for the body-beam-body 
configuration for "pitch" bending, "roll" bending and "yaw" twisting. 
Since the modes are based on solving explicitly the distributed parameter 
equations (without damping and without kinematic coupling) there is no 
limit to the number of modal characteristic sets that can be generated by 
the program. It will be the analyst's decision as to how many modes need 
to be considered. 


Laboratory Experiment Description 

The second part of the design challenge is to validate in the 

laboratory, the system performance of the more promising control system 
designs of the first part. The experimental apparatus will consist of a 
dynamic model of the Space Shuttle orbiter with a large antenna reflector 
attached by means of a flexible beam. The dynamic model will be exten- 
sively instrumented and will have attached force and moment generating 
devices for control and for disturbance generation. A single, flexible 
tether will be used to suspend the dynamic model, allowing complete angular 
freedom in yaw, and limited freedom in pitch and roll. An inverted 

position will be used to let the reflector mast to hang so that gravity 

effects on mast bending will be minimized. The dynamics of the laboratory 
model will of necessity be different from the mathematical model discussed 
earlier. 
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Design Challenge, Part One 


For part one of the design challenge, the following mathematical 
problem is addressed. Given the dynamic equations of the Shuttle/antenna 
configuration, what control policy minimizes the time to slew to a target 
and to stabilize so that the line-of-sight (LOS) error is held, for a time, 
within a specified amount, <$. During the time that the LOS error is 
within 6, the attitude must change 90° to prepare for the next slew 
maneuver. This was previously referred to as the sescondary control task. 
The maximum moment and force generating capability will be limited. Advan- 
tage may be taken of selecting the most suitable initial alignment of the 
Shuttle/antenna about its assigned initial RF axis, line-of-sight. 

Random, broad band-pass disturbances will be applied to the configuration. 
Two proof-mass, force actuators may be positioned anywhere along the beam. 
The design guidelines are summarized below: 

1. The initial line-of-sight error is 20 degrees. 

e(o) = 20 degrees 

2. The initial target direction is straight down. 



3. The initial alignment about the line-of-sight is free to be chosen 
by the designer. Advantage may be taken of the low value of 
moment of inertia in roll. The Shuttle/antenna is at rest 
initially. 

A. The objective is to point the line-of-sight of the antenna and 
stabilize to within 0.02 degree of the target as quickly as 
possible. 

5 * 0.02 degree 
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5. Control moments can be applied at 100 Hz sampling rate to both the 
Shuttle and reflector bodies of 10,000 ft-lb for each axis. The 

commanded moment for each axis is limited to 10,000 ft-lb. The 

actual control moment's response to the commanded value is 
first-order with a time constant of 0.1 second. 

For the rolling moment applied to the Shuttle body: 


-10 < 


- Yl, 


. < 10 
command — 


- 0.1 


- 0 . 1 , 


Yl*” + 0 - e ' Yl 00 + (1 - * > Yl,co«..nd<"> 


Equations for other axes and for the reflector body are similar. 

6. Control forces can be applied at the center of the reflector in 
the X and Y directions only. The commanded force in a 

particular direction is limited to 800 lbs. The actual control 
force's response to the commanded value is first-order with a 
response time of 0.1 second. 

For the side for applied to the reflector body: 


-800 < F . < 800 

— Y, command — 


F y (n + 
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e" 0 * 1 F y (n) 


+ (1 " e" 0,1 ) 


F v .(n) 

Y, command 


Equations for X-axis are similar. 

7. Control forces using two proof-mass actuators (each having both 
X and Y axes) can be applied at two points on the beam. The 

strokes are limited to ± 1 ft, and the masses weight 10 lbs each. 
The actual stroke follows a first-order response to limited 


2 SO 


commanded values* 


For the X-axis of the proof-mass actuator at s 2 : 


-1 < A . < 1 

— X ,2, command — 


4 X,2<" + »> * e '° -1 4 X.2<”> + **.*.«—« 


Equations for other axes and locations are similar. 

8. The inertial attitude direciton cosine matrix for the Shuttle body 
lags in time the actual values by 0.01 second and are made at a 
rate of 100 samples per second. Each element of the direction 
cosine measurement matarix is contaminated by additive, 
uncorrelated Gaussian noise having an rms value of 0.001. The 
noise has zero mean. 



d l 1 ( n ) d 12 ( n ) 

d 13 

T ,(n + 1) = T _ (n) + 

s, measured s,true 

d 21^ d 22^ 

d 23 


d 31 (n) d 32 (n) 

d 33 

where : 



Ejd^Cn)} = 0 

E{d 1J (n)d kL (n)} - 0 
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for k * 0 
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9. The angular velocity measurements for both the Shuttle and 
reflector bodies pass through a first-order filter with 0.05 sec 
time constant and lag in time the actual values by 0.01 second and 
are made at a rate of 100 samples per second. Each rate 

measurement is contaminated by additive, Gaussian, uncorrelated 
noise having an rms value of 0.02 degree per second. The noise 
has zero mean. 

For example: 


(0 (n + 1) a u) 

1 ,X, measured' 1,X, filtered 

E { e i x (n) Cj x (n + k)} - 0 

• (. 02) 2 


(°> + e ljX („> 
for k * 0 

for k * 0 


where 


to = - 20 to + 20 to 

1,X, filtered 1,X, filtered l,X,true 


10. Three-axis accelerometers are located on the Shuttle body at the 
base of the mast and on the reflector body at its center. Two- 
axes (X and Y) accelerometers are located at intervals of 
10 feet along the mast. The acceleration measurements pass 
through a first-order filter with a 0.05 second time constant and 
lag in time the actual values by 0.01 second, and are made at a 
rate of 100 samples per second. Each measurement is contaminated 
by Gaussian additive, uncorrelated noise having an rms value of 
0.05 ft/sec^. 
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For example: 


a l,X,measured (n + 0 “ a i ,x,filtered (n) + T l,X (n) 

E ^ T 1 X (n) T 1 X (n + k >l “ 0 for k * 0 

- (.05) 2 for k = 0 

where : 

^1,X, filtered " 20 3 1 ,X, filtered + 20 “l.X.true 

11. Gaussian, uncorrelated step-like disturbances are applied 
100 times per second to the Shuttle body in the form of 3-axes 
moments, having rms values of 100 ft-lbs. These disturbances 
have zero mean. 

For example: 

E {M D x (n) M D>x (n + k)} = 0 for k ^ 0 

= ( 100) 2 for k = 0 

In summary, the designer’s task for part one is to: (1) derive a 

control law for slewing and stabilization, coded in FORTRAN; (2) select an 
initial attitude in preparation for slewing 20 degrees; and (3) select two 
positions for the 2-axes proof-mass actuators. An official system 
performance assessment computer program will be used to establish the time 
required to slew and stabilize the Shuttle/antenna configuration. 


253 


Design Challenge, Part Ttoo 


As in part one, the task is to minimize the time to slew and stabilize 
a Shuttle/antenna configuration. The difference is that in part two of the 
design challenge, a physical laboratory model will be used instead of the 
dynamic equations of part one. The constraints on total moment and force 
generation capability will apply to part two, as for part one. Again, the 
analyst may select the initial alignment about the assigned initial RF 
line-of-sight. Disturbances will be injected into the Shuttle/antenna 
model. The designer's task will be similar to that for part one. 

CONCLUDING REMARKS 

A Design Challenge, in two parts, has been offered for the purpose of 
comparing directly different approach to controlling a flexible 
Shuttle/antenna configuration. The first part of the design challenge uses 
only mathematical equations of the vehicle dynamics; the second part uses a 
physical laboratory model of the same configuration. The Spacecraft 
Control Laboratory Experiment ( SCOLE) program is being conducted under the 
cognizance of the Spacecraft Control Branch at the NASA Langley Research 
Center. The NASA/ IEEE Design Challenge has the advice and counsel of the 
IEEE-COLSS Subcommittee on Large Space Structures. Workshops will be held 
to enable investigators to compare results of their research. 
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Figure 1. Drawing of the Shut tie /Antenna Configuration. 
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Figure 2.- Schematic of the effect of bending on the 
line-of-sight pointing error. 
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CONTROL MOMENTS, DISTURBANCES 




wunoe AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE AMPLITUDE 


Roll bending mode # 1. Frequency= .32 Hz 



Roll bending mode # 2. Frequency= 1.29 Hz 



Roll bending mode § 3. Frequency= 480 Hz 



Roll bending mode § 4 Frequency= 12.29 Hz 



Roll bending mode # 5. Frequency= 23.68 Hz 



Roll bending mode § 7. Frequency = 57.90 Hz 



Roll bending mode § 8. Frequency= 80.72 Hz 
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Figure An .- Plots of normalized roll bending mode shapes 
for SCOLE configuration. 
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Pitch bending mode # 1. Frequency= .29 Hz 



Pitch bending mode # 2. Frequency= 1.65 Hz 



Pitch bending mode § 3. Frequency= 4.97 Hz 
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Pitch bending mode # 4 . Frequency= 12.36 Hz 



•«F 




Pitch bending mode 0 7. Frequency 57.92 Hz 

P 

Pitch bending mode § 8. Frequency= 80.73 Hz 



Figure 4b*— Plots of normalized pitch bending mode shapes 
for SCOLE configuration. 
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Torsional mode # 1. Frequency = .53 Hr 



Torsional mode # 2. Frequency = 45.12 Hr 



Torsional mode # 3. Frequency= 90.23 Hz 



Torsional mode $ 4. Frequency— 135.35 Hz 



Torsional mode # 5. Frequency= 180.46 Hz 



Torsional mode # 6. Frequency= 225.57 Hz 



Torsional mode # 7. Frequency= 270.69 Hz 



Torsional mode # 8. Frequency= 315.80 Hz 


Figure Ac.- Plots of normalized torsional mode shapes for 
SCOLE configuration. 
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